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Abstract 

An extension of the generalized relativistic mean-field (gRMF) model with 
density dependent couplings is introduced in order to describe thermodynam- 
ical properties and the composition of dense nuclear matter for astrophysical 
applications. Bound states of light nuclei and two-nucleon scattering corre- 
lations are considered as explicit degrees of freedom in the thermodynamical 
potential. They are represented by quasiparticles with medium-dependent 
properties. The model describes the correct low-density limit given by the 
virial equation of state (VEoS) and reproduces RMF results around nuclear 
saturation density where clusters are dissolved. A comparison between the 
fugacity expansions of the VEoS and the gRMF model provides consistency 
relations between the quasiparticles properties, the nucleon-nucleon scatter- 
ing phase shifts and the meson-nucleon couplings of the gRMF model at zero 
density. Relativistic effects are found to be important at temperatures that 
are typical in astrophysical applications. Neutron matter and symmetric 
nuclear matter are studied in detail. 
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1. Introduction 



Correlations are an essential feature of an interacting many-body system 
such as nuclear matter. They have an impact on the thermodynamical prop- 
erties and the composition of the matter since the constituent particles can 
form new nuclear species as bound states. This aspect is particularly im- 
portant in the description of matter below nuclear saturation density that is 
employed in theoretical simulations of astrophysical phenomena such as core- 
collapse supernovae. There, the knowledge of the equation of state (EoS) of 
neutron-proton asymmetric matter is required in a wide range of densities 
and temperatures [H, The description of dense matter is important for 
the investigation of various stages in supernova explosions and the structure 
of neutron stars [3]. It is known to affect the effectiveness of the neutrino 
reheating of the shock wave [4]. The sensitivity of the collapse dynamics on 
the properties of matter and the composition in this density regime could 
strongly influence the structure of the proto-neutron star 

Although there are many approaches to describe dense matter, in par- 
ticular at zero temperature, only few models cover the full parameter space 
needed in most astrophysical applications. Often, they do not supply suffi- 
cient information on the thermodynamical and compositional details. Thus, 
for many years, a very small number of EoS tables was available which have 
been used in simulations of dynamical astrophysical processes How- 
ever, in recent years, the interest in developing EoS has surged and several 
new EoS tables were provided |9-16|. This development was triggered by the 
needs of astrophysical modelers, the progress in the theoretical description 
of nuclear matter and an increase of computational power. Nevertheless, 
approximations and simplifications are still needed for practical purposes. 

There are two major paths to build an EoS of warm and dense nuclear 
matter for practical applications. One approach starts with an ideal mix- 
ture of nucleons and nuclei leading to a nuclear statistical equilibrium (NSE) 
description \v\\. The effect of interactions between all constituents can be 
incorporated with the help of viriai corrections. However, at present, only 
nucleons and light nuclei are considered in practice in such a viriai equation of 
state (VEoS) |l7Hl9l| that provides the correct finite-temperature EoS in the 
limit of very low densities. The results are model independent since they de- 
pend only on experimentally determined data, i.e. binding energies of nuclei 
and scattering phase shifts. Unfortunately, the application of this approach is 
limited to rather low densities. The dissolution of nuclei and the transition to 
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uniform neutron-proton matter with increasing density cannot be described 
properly. In order to simulate such an effect, the heuristic excluded vol- 



ume mechanism was employed frequently A second class of EoS models 



for astrophysical applications is based on self-consistent mean-field methods 
with neutrons and protons as fundamental constituents. They are consid- 
ered as quasiparticles with self-energies that contain the information on the 
interaction which is usually modeled in an effective way and not taken from 
a realistic nucleon-nucleon (NN) interaction. These mean-field models, both 
non-relativistic and relativistic, can be very successful in describing finite 
nuclei and nuclear matter around saturation density rzsat ^ 0.16 fm 



But at low densities n <^ n^ax. they fail to include correlations properly that 
give rise to inhomogeneities and the formation of many-body bound states, 
i.e. nuclei. There were attempts to modify the low-density behavior of mean- 



field models guided by microscopic calculations, see, e.g., Ref. [2J] for the 
zero temperature case. However, it is still challenging to include the forma- 
tion of clusters properly. The deficiencies of the various models lead to the 
strategy of patching up or merging different approaches providing a uniform 
description of warm dense matter. 

A recent publication ^] considered a quantum statistical approach to nu- 
clear matter and devised a generalized relativistic mean-field (gRMF) model, 
which is an extension of the relativistic mean-field model with density de- 



pendent couplings (DD-RMF) from Ref. [26||. The gRMF model includes, 
besides nucleons, light nuclei (^H = d, = t, ^He = /i, ^He = a) as hadronic 
degrees of freedom in the Lagrangian. The dissolution of these clusters is 
modeled by a medium-dependent shift of their binding energies originating 
mainly from the action of the Pauli principle. All hadronic constituents of 
the model can be viewed as quasiparticles that interact via the exchange 
of effective mesons. The gRMF model achieves bridging the two views on 
nuclear matter, i.e. a mixture of nucleon and nuclei at low densities and nu- 
cleons as quasiparticles moving in mean fields at high densities. Even though 
the model contains the same relevant particles at low densities as the VEoS, 
it does not reproduce exactly the thermodynamical properties of the VEoS 
as will be shown below. In fact, all presently employed EoS for astrophys- 
ical applications that are based on mean-field concepts share the problem 
of reproducing the VEoS in the low-density limit. Thus the question arises, 
how the models can be modified in order to approach the correct low-density 
limit. Is a modification of the effective interaction sufficient to meet this aim 
or are there more severe changes of the models required? 
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In this work we propose an approach that introduces bound and scatter- 
ing states in an effective way as additional degrees of freedom in the ther- 
modynamical potential. All relevant quantities are derived in a thermody- 
namically consistent way. These many-body bound and scattering states are 
represented by quasiparticles with medium-dependent properties and tem- 
perature dependent energies. The low-density behavior of nuclear matter 
at finite temperatures is considered in detail by comparing the VEoS with 
the gRMF approach by means of a series expansion of the grand canonical 
potential in powers of the nucleon fugacities. From the comparison of the 
VEoS and gRMF expansions we will derive consistency relations that con- 
nect quasiparticle properties with the meson-nucleon couplings in the vac- 
uum and the phase shifts or effective-range parameters of nucleon-nucleon 
scattering. We investigate different choices of the meson-nucleon couplings 
and quasiparticle properties set by these consistency relations. In addition, 
relativistic corrections to the traditional VEoS are obtained from these re- 
lations that become larger than the effects of particle correlations in the 
low-density limit. Various parametrizations of the quasiparticle medium de- 
pendence are discussed. Studying thermodynamical properties of the matter 
we consider temperatures higher than the critical temperatures for pairing 
and Bose condensation in order to avoid these more complex effects. Thus 
our consideration of the low and zero temperature limits should be corrected 
in the future to incorporate condensation phenomena. Our approach is not 
limited to the chosen particular gRMF model with density-dependent meson- 
nucleon couplings. It can be applied easily to other mean-field approaches 
as well. 

This paper is organized as follows. In Section [2] we present a general 
formulation of the VEoS that serves as one starting point of the discus- 
sion. Basic quantities are defined and the notation is established. We show 
that correlations in continuum states can be represented effectively by tem- 
perature dependent resonances. Analytic expressions for the second virial 
coefficients are found with the help of the effective-range expansion for the 
s-wave phase shifts. The connection of the VEoS and NSE models is estab- 
lished. The relation between the VEoS and the generalized Beth-Uhlenbeck 
approach, which takes into account modifications due to medium effects, is 
discussed. The gRMF model with density-dependent meson-nucleon cou- 
plings is sketched in Section [3] and a series expansion of the grand canonical 
potential Q in powers of fugacities is derived for low densities. The power 
series for Q in the VEoS and gRMF are compared up to second-order leading 
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to consistency relations between the models that are studied in various limits. 
In Sect. IHwe consider an extension of the gRMF model: in order to satisfy 
consistency relations we introduce temperature dependent resonance energies 
of the scattering states and modified degeneracy factors that allow to make 
a smooth interpolation between VEoS and gRMF models. We find that rela- 
tivistic corrections to the standard VEoS are important in these consistency 
relations already in the first order of the expansion. The particular example 
of neutron matter is presented in Section [5] and symmetric nuclear matter is 
discussed in Section [61 The following section deals with the transition from 
low to high densities and the occurring problems. Concluding remarks and 
an outlook are given in Section |H1 The comparison of our notations with the 



ones given in Ref. [18| can be found in Appendix A An expansion of the 



energy per particle in neutron matter at zero temperature in powers of the 



Fermi momentum is considered in Appendix B Throughout this work we 
use natural units such that h = c = 1. 

2. Equation of state in the virial limit 

2. 1 . General formalism 

The VEoS represents a model-independent approach in the calculation of 
thermodynamical properties of low-density matter at finite temperatures T. 
In the non-relativistic limit, which is usually considered, it describes a system 
of interacting particles i, j, . . . with non-relativistic chemical potentials in 
a volume V , provided the fugacities Zi = exp{fii/T) are small {zi ^ 1). The 
range of densities where the virial expansion is valid can be estimated by the 
relation njAf <C 1 where is the number density of particle i with mass m,i 
and Aj = ^/27r/ {miT) is the thermal wavelength. 

The following presentation gives a generalized and more symmetric for- 



mulation of the virial approach as compared to [17h19|. The S- matrix for 



mulation by Ref. [27j that was recently applied to nuclear matter in Ref. [28 
gives essentially identical results. 

Under the presupposed conditions the grand canonical partition function 
Q can be expanded in powers of fugacities as 

Q(T, V, /ii) = 1 + X QiZi + ^ QijZiZj + ^ XI QijkZiZjZk + ... (1) 
i ij ijk 

with one-, two-, three-, . . . many-body canonical partition functions Qj, Qjj, 
Qjjfc, .... In classical non-relativistic mechanics, the single-particle canonical 
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partition function is 

Qi-J^sJ d'ri j Spi exp {-(iH,) = (2) 

witli tlie spin degeneracy factor gi (= 2 for each nucleon), (3 = 1/T and the 
single-particle Hamiltonian — pj/ {2mi) that only contains a kinetic con- 
tribution with momentum pi. For the two-body canonical partition function 
we have 

Q^j = J d'n J d^'p, j j d'p, cxp{-/3H,,) (3) 

with the two-body Hamiltonian Hij — pf/ {2mi) +Pj / {2mj)+Vij that includes 
a two-body potential Vij responsible for the correlations. 
Then, the grand canonical potential 

Q(T, V, lii) = -T In Q(T, 1/, = -pV , (4) 

that is directly related to the pressure p, can be written in the form 

n{T,V,f,,)= (5) 

with the (symmetrized) dimensionless cluster (or virial) coefficients 



bi{T) = ^g,, (6) 

^3/2^3/2 

bij{T) = ^_^(Q.._Q.g.) , (7) 

bijk{T) = ^ {Qijk — QiQjk — QjQik — QkQij + ^QiQjQk) (8) 

that depend on the temperature. Without interaction, the two-, three-, 
...many-body partition functions factorize, e.g. Qij — QiQj, and the sec- 
ond, third, . . . cluster coefficients bij, bijk, ■ ■ ■ vanish. 

Individual particle number densities are found from the relation 



1 dQ 
V di^i 
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with contributions from free particles (first term in the second equation) and 
correlated pairs, triples, etc. The entropy is 



S 



on 

df 

5n 



(10) 



with coefficients 



A? 



A 



•A- 



ijk 



XiXjXk 



dT T ' 



^ijk 



T- 



dh 



T 



dT T ' 



T 



(11) 
(12) 

(13) 



Other relevant thermodynamical quantities such as the free energy and the 
internal energy 

F = n + vY,^^^^i (14) 



and the internal energy 



E = F + TS 



(15) 



can be obtained immediately. 

The formulas given above can be generalized to include four-, five-, . . . many- 
body correlations, but already contributions of the three-body term are 
hardly ever considered in practice. For a given temperature, two-body cor- 
relations will always dominate higher-order correlations when the density 
decreases. Hence, in the following, we will truncate the expansion at second 
order in the fugacities of the basic constituents. 

The quantum mechanical generalization of the virial expansion up to sec- 
ond order was given by Beth and Uhlenbeck 29|, l30| . In classical mechanics 
the interaction potential between the particles is the relevant quantity that 
appears in the calculation of the virial coefficients. In quantum mechanics, 
the Schrodinger equation has to be solved with this potential and the ob- 
tained eigenstates comprise bound and scattering states. Integrations over 
phase space are replaced by sums over all eigenstates of the one- and two- 
body system and it is obvious that the density of states becomes the relevant 
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quantity. The result for the first virial coefficient bi is identical to the classical 
value in the continuum approximation, corresponding to the replacement 
of the sum over discrete momentum states by an integral. The second virial 
coefficient can be expressed as an integral over center-of-mass energies E 

hAT) = \ ' ^3,' j dE exp i-^E) D,,{E) ± 6,^9,2-'^' (16) 
with \ij = ^J2tx j\{mi + vn^)T\ and the quantity 

A. (£;) = ^?'^(^ - ^) + E (17) 

k I 

that measures the difference between the density of states of an interacting 
and free two-particle system. The last term in equation ffT6l) is a quantum 
statistical correction with the positive (negative) sign if i = j are identical 
bosons (fermions). 

The first contribution to is a sum over all two-body bound states k of 
the system ij with degeneracy factors gi.^'^ and energies Ef'^ < 0. The second 
term takes into account continuum correlations in the two-particle system in 
all channels / with degeneracy factors g^^'^^ via the energy-dependent scatter- 
ing phase shifts ^f"'''. In general, the index / stands for the total and angular 
orbital momenta of a particular channel. If the phase shift is dominated by 
a narrow resonance, it jumps by vr in a very small energy interval around 
the resonance energy and the contribution to the virial coefficient resembles 
that of a bound state. With experimentally known bound state energies and 
phase shifts, the second virial coefficient and thus the low-density EoS 
can be calculated in a model-independent way. 

The quantum mechanical generalization of the third virial coefficient was 



discussed in Ref. [31| but actual calculations turn out to be difficult in prac- 
tice and, thus, they are not considered here. The derivation of the classical 
VEoS is based on non-relativistic kinematics. In order to include relativistic 
effects, at least the relativistic dispersion relation of the particles has to be 
used. The resulting modifications of the virial coefficients will be discussed 
when the VEoS is compared to the gRMF approach in Section 13.41 

2.2. Application to nuclear matter with arbitrary neutron to proton ratio 

In low-density nuclear matter at not too high temperatures, neutrons 
and protons are the relevant fundamental constituents. Hence, the grand 
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canonical potential up to second order in the fugacities becomes 



(18) 



rp-rr I ~n , "p , "nn 2 I "pp 2 
I A3 ^ A3^^ ^ ^""^ ^ A3 



''np 



' ,3/2, 3/2 ^"^P 
An Ap 



where the symmetry = 6„„ is used. The Coulomb interaction is usually 



neglected in the VEoS and 6„„ is replaced by bnn in ( 1181) . cf. Ref. [18[. The 



second virial coefficient bnn receives contributions only from scattering states. 
It is convenient to split bnp = &npo + ^npi into two contributions with total 
isospin T = and T = 1. The only bound two-body state in the two-nucleon 
system, the deuteron, appears in the isospin zero channel with (7^ = 3 and 



E. 



{np) 



bnn{T) 

bpp{T) 

bnpl{T) 
bnpO (T) = 

with the virial integral 



2.225 MeV [32(]. According to this observation we can write 



A^ 



nn 
3 



P \^ App) t{pp) _ ^ 0-5/2 
\3 Z^^l 



A3 
PP 



1A: 
2 



3/2,3/2 
Ap 



A3 

np 



E (npl) ^(npl) 



3/2^3/2 
P 



ia;^'^a 

2 



A3 

'np 



S-d exp 



(npO) ^(npO) 



irm = — 



dE d6. 



I 



TT 



dE 



exp — — 



T 



Formally, we can write the sum of the scattering contributions in Eqs. 
fl22|) in the form of a single bound state contribution 



(19) 
(20) 
(21) 
(22) 

(23) 

m- 



°° dE d6^j 



TC 



dE 



exp 



E 
T 



9ij exp 



E^J{T) 

T 



(24) 



with a temperature dependent effective resonance energy Eij{T) in each 
channel by summing the contributions of all partial waves in an effective total 
phase shift 

(25) 



S.,{E) 
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The effective degeneracy factor cjij = ifi'o*"''' is chosen such that it is identical 
to the degeneracy factor in the s-wave. The sign is determined by that of the 
integral. 

In order to illustrate the behavior of the phase shifts and effective reso- 
nance energies, neutron-neutron scattering is considered in the following. In 
this example, the sum of the experimental phase shifts receives contributions 
only from certain partial waves with total angular momentum J and orbital 
angular momenta / due to the Pauli principle. Taking into account all pos- 
sible partial waves with / = 0, 1 and 2, i.e. the channels ^5*0, ^Pq, ^Pi, ^P2 
and in spectroscopic notation ^^^^Lj, one has 

6^n{E) = 5^7^ + 5(7) + 3(5ir^ + 5(5g") + S^g") (26) 

using double indices Jl in 6jp to indicate the partial wave. In Fig. [1] the 
effective phase shift for the sum of all partial waves with / < 2 is depicted 
in comparison with the pure s-wave phase shift. Since there is no exper- 
imental data on the nn phase shifts, we use the np data taken from the 
Nijmegen partial wave analysis [33| to illustrate the comparison between the 
experimentally known phase shifts and the results using the effective range 
approximation. Note, that the laboratory energy i?iab = 2i? of a neutron 
scattered on the target at rest is used as the argument. At very low energies, 
both curves are rapidly rising. At energies above ~ 10 MeV, contributions 
of higher partial waves become important. 

At low temperatures, the virial integral is dominated by the low-energy 
phase shifts 61 and only the s-wave contributes significantly due to the large 



derivative with respect to the energy. The effective-range approximation [34 
for the s-wave phase shift 

kcot5l;'^ = -— + lnje (27) 

CLij / 

with momentum k = ^j2iJLijE and effective mass /ijj = mimj/{mi + m-,) can 

be used to specify the energy dependence of "'^ with help of the scattering 
length aij and the effective range r^j. 

Fig. [1] shows the result for the s-wave phase shift using only the scattering 
length or both the scattering length and the effective range parameter. The 
experimental phase shift is nicely reproduced by the effective range expansion 
at very low energies. The contribution of the fc^ term with rij is essential to 
obtain a decrease of the s-wave phase shift at higher energies. In this case. 
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Figure 1: Effective phase shift (5„„(£'iab) for neutron-neutron scattering versus laboratory 
energy i?iab = '2E. Resuhs are shown using experimental data (including s, p and d-waves 
or s-waves only) and the effective-range approximation with and without the contribution 
depending on the effective range parameter, respectively. 



the result follows more closely the experimental data, however, differences 
with respect to the total effective phase shift 6nn{E) remain. 

The effective resonance energy Enn for nn scattering derived from the 
phase shifts in Eq. ( 12^ is compared for different approximations in Fig. |2j 
In general, Enn rises smoothly with increasing temperature indicating that it 
is not dominated by a particular resonance. Contributions of higher partial 
waves reduce the effective resonance energy as compared to the pure s-wave 
result. 

The s-wave virial integral can be obtained analytically in the effective- 
range approximation as 



5, 



iv) 



exp 



for 2rij/aij < 1. It depends on the coefficients 

,2 



A. 



iv) 



2 



erfc 



1 



(28) 



(29) 
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Figure 2: Temperature dependence of the effective resonance energy i5„„ for neutron- 
neutron scattering in different approximations. See caption of Figure [Hand text for details. 



and 



l + r/./l-2^ 



For = the integral reduces to 



(30) 



2 a,; 



exp 



2fiijTa: 



erfc 



(31) 



Using these results, the effective resonance energy can be expressed explicitly 
as a function of the temperature and the effective-range parameters. We 



use ann = —18.818 fm and r„„ = 2.834 fm from Ref. [3a] . As shown in 
Figure [2], Enn in the effective range approximation rises more strongly with 
temperature when the contribution of the r„„ term is taken into account. 
At high temperatures there remains a difference to the exact result that 
takes higher partial waves into account. However, by a readjustment of the 
scattering length to a„„ = —11.31 fm and the effective range to r„„ = 1.06 fm, 
the correct effective resonance energy obtained with all partial waves is very 
well reproduced for the relevant temperatures with the modified parameters 
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(denoted by modified a„„ + r„„ in Figure [2]). Of course, tliese modified 
values for a„„ and r„„ are not pliysical and do not describe tlie scattering in 
a particular partial wave. Tfiey are only an effective means to describe the 
temperature dependence of the effective resonance energy. 

There are several interesting limits for the virial integral. For low tem- 
peratures T ^ A[^^ /{2ij,ij) one finds 



4-)(r)^-a,Y^ + ... (32) 

(independent of the effective range Vij). In the unitary limit aij — oo 
and = 0, that describes the situation where a bound state/resonance is 
just at the continuum threshold, Iq^\t) = 1/2 and the effective resonance 
energy is given by Eij = Tin 2. This result is also depicted in Fig. [2l For 
positive scattering length aij, e.g. in the deuteron channel, the integral is 
negative and compensates part of the correlation strength that is located in 
the bound state of the same channel. 



In Ref. [18| also the a-particle was considered as a fundamental con- 



stituent in the VEoS and the additional virial coefficients 6q,, 6„q,, bpa and 
baa appeared in the formalism. (For a comparison of the virial coefficients 



in the two formulations see Appendix A ) In our approach, the correspond- 



ing contributions to represent four-, five- and eight-body correlations that 
will become irrelevant in the low-density limit as compared to two-nucleon 
correlations. Since chemical equilibrium requires jj^ = 2/i„ -|- 2fip + Ba with 



the a-particle binding energy Ba = 28.296 MeV [32|, the a-particle fugacity 
Za = exp{fia/T) is not independent of the neutron and proton fugacities. The 
VEoS with neutrons, protons and a-particles was further extended in Ref. 
[l9] to include and ^He as constituent particles corresponding to three- 
body correlations. Care has to be taken in order to avoid a double counting 
of states since, e.g., the a-particle states can be obtained in the ^H-p and 
^He-n channels (or the neglected ^H-^H channels) that all represent four-body 
correlations when nucleons are assumed to be the only fundamental particles. 

2.3. Relation to the nuclear statistical equilibrium approach 

In principle, the virial expansion could be extended to third, fourth, 
. . . order to include more and more many-body correlations. The second, 
third, . . . virial coefficients will contain bound state contributions that cor- 
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respond to ground and excited states of heavy nuclei. Neglecting the many- 
body continuum correlations, the grand canonical potential becomes a sum 



s^^-rv-Et^-pf^) (33) 



over all nuclei with mass number A, charge number Z, chemical potential 
= (A — Z)fj,n + ZfXp and thermal wavelength Xa,z = A/27r/ {■mA,zT). 
The effective degeneracy factor 

9aAT) = 9% + E 9^AZ (-f ) (34) 

is a sum over the ground state and all excited states k of the nucleus 
with excitation energies E^, total angular momenta Jk and g^^\ = {2Jk + 1). 
The summation over individual states is often replaced by an integral over 
energy with an appropriate level density, see, e.g. Ref. The temperature 
dependence of the degeneracy factor is a simple and efficient means to include 
excited states of nuclei. It will be used below in the formulation of the gRMF 
approach. 

2.4- Generalized Beth-Uhlenbeck approach 

The description of matter in the conventional VEoS is limited to rather 
low densities. However, the approach can be generalized in order to reach 
higher densities and to include the effects of dissolution of composite parti- 
cles, i.e. the Mott effect. This generalized Beth-Uhlenbeck (gBU) approach 
is based on a quantum statistical description with thermodynamic Greens 



functions. For details see Ref. [36[. The grand canonical potential f2 in 
the gBU approach assumes a similar form as in the VEoS. The constituents 
are considered as quasiparticles with the correct statistics (Fermi-Dirac or 
Bose-Einstein) and with self-energies that contain already some effect of the 
mutual interaction. The remaining two-body correlations (between quasi- 
particles!) are contained in a modified second virial coefficient that is given 
by 

KiT) = ^-^Xfxf j dE U,{E + E,,^,)F,,{E) (35) 

with the correct distribution function fij{E) of the composite system, i.e. 
Bose-Einstein for the two-nucleon-states, and the quantity 
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F.AE) = Y.^' I M^J^E-^t'') (36) 
+ 2^ ^ J (2„)' ^ ' dE 



that is related to the in-medium density of states. The main differences as 
compared to Eqs. ( ITB]) and (fT7|) are caused by medium effects. The prop- 
erties of a two-body system ij depend on the total cm. momentum P of 
the state with respect to the medium. Hence, there is an additional sum- 
mation (integration) over P. Bound states only appear for P larger than 
the Mott momentum Piviott because at lower total momenta their forma- 
tion is suppressed by the action of the Pauh principle. The phase shifts 
fii,fij) are determined by the in-medium T-matrix describing the 
scattering of two quasiparticles. Similarly Ej^-'\P,T, jj^i, jj^j) is the medium- 
dependent energy of a bound state. Energies in Eq. fl36|) are measured with 
respect to the continuum edge -Econt that corresponds to the energy in the 
frame of the medium where the relative momentum of the scattering parti- 
cles is zero. The 2 sin^ 6^^''^ factor takes into account that part of interaction 
effects is already contained in the self-energies of the quasiparticles. The 
second term in Eq. (!36|) is related to the number of correlated quasiparticles 
that is different to the number of correlated particles as calculated by the 
usual Beth-Uhlenbeck approach in Eq. ([1] 



3. Generalized relativistic mean-field model 

3. 1 . General formalism 

In Ref. M a RMF model with density-dependent meson-nucleon cou- 



plings was generalized to include light clusters (bound states of ^H, ^H, ^He, 
^He) as additional degrees of freedom with medium-dependent binding ener- 
gies. These clusters couple minimally to the meson fields with a strength that 
was assumed to be a multiple of the nucleon coupling strength. Scattering 



states were not considered in the RMF model of in Ref. [25|]. The approach 
can easily be extended when more than the above mentioned nuclei need to 
be included. Such a model with nuclei as explicit degrees of freedom we call 
the gRMF model in this paper. 
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The grand canonical potential of the model can be formulated as a func- 
tion fl(T,V, p,a,6) depending on the temperature T, volume V, rela- 
tivistic chemical potentials fii = rrii + fii of the particles and the spatially 
constant meson fields co, p, a and 6 under consideration. Here, all four mesons 
are included to cover the four possible combinations of scalars and vectors in 
Lorentz and isospin space. We neglect antiparticles in the following because 
they give sizable contributions only at temperatures much higher than those 
relevant in the considered astrophysical applications. In the general case of 
inhomogeneous matter, not discussed here, fl becomes a functional of the 
spatially varying meson fields and their gradients. 

For uniform nuclear matter, the grand canonical potential can be written 



as 



n 



V 



^ / 2 2 , 2 2 2 2 



(37) 



+ {T'^ujn^ + T'ppnp - T'^an^ - T's5ns) (n„ + 



where F- are the derivatives of the density dependent couplings which de- 
termine the nucleon-meson coupling strength. Contributions from individual 
particles (nucleons and clusters) are 



ili = TQiVT 



(2^ 



In 



1 ± exp 



Ej — Pi 
T 



(38) 



where the upper (lower) signs refer to fermions (bosons). In the case of 
Maxwell-Boltzmann statistics, Eq. fl38|) reduces to 



= g,VT 



(2^ 



exp 



T 



(39) 



We assume that the degeneracy factors gi can depend in general on the 
temperature as in the case of NSE models, cf. Section [2731 The densities n^. 



n 



pi 



Her and ns appearing in ( 137|) are themselves functions of the temperature, 
chemical potentials and meson fields (see below). 

The mass of a cluster i = d,t,h,a, . . . with Ni neutrons and Zi protons 
is written as 



rrii 



NilTLn + ZiTTlp - (1 - 6in){l - ^ip)Bl 



(vac) 



(40) 



with neutron and proton rest masses m„ and rrip, respectively, and the vac- 
uum binding energy > for composite particles. For neutrons and 
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protons the binding energy does not appear. Chemical equihbrium leads to 
the constraint 

jii = Nijln + ZifXp (41) 

for the relativistic chemical potentials of the clusters. With every particle 
i = n,p, d, t,h,a, . . . of momentum k, an energy 

Ei{k) = V,+ + (m, - S,y (42) 

is associated that contains the vector and scalar self-energies Vi and Si, re- 
spectively. These potentials are given by 

V, = Ti^io + T,,p+{6,n + S^p)V^'■\ (43) 
Si = Ti^a + T,s6-il-6,n){l-5^p)ABi (44) 

with meson fields m = co, p,a,6 and their couplings Tim to the particle i. We 
assume that the coupling strength of the clusters is a multiple of that of the 
nucleons, i.e. 

^im = girrJ^miQ) (45) 

with 

giuj = 9ia = Ni + Zi, (46) 
gip = QiS = Ni- Zi. (47) 

Other choices are possible and should be explored in the future. The meson- 
nucleon couplings Tm{Q) depend on the total nucleon density g = Un + rip 
only. The functional form and the parameters of the density dependent 
couplings Tmig) used in the gRMF model were chosen in such a way to 
describe the properties of finite nuclei and nuclear matter parameters at 
saturation density. The vector self-energy contains the rearrangement term 

V^^*-) = Tin^u + T'^n^p - T'^n^a - T',ns5 (48) 

that is only relevant for nucleons. In the medium, the actual binding energy 
of a cluster Bi = B^^^^'' — ABi will be shifted with respect to the vacuum 
value B^^^^^ by a quantity ABi ■ This medium-dependent binding energy shift 
ABi for composite particles is parametrized as a function of temperature and 
the Lorentz vector meson fields u and a, see, e.g. Ref. [25j. The specific form 
of these shifts is not important in the present discussion of the low-density 
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behavior of the EoS but it affects the transition to high densities, see Section 

m 

The source densities in Eqs. f l37|) and fHHj) 

= ^ Qiurii , np = ^ QipTii , (49) 

(50) 



giaU] ' , ns = 2_^ gisnl 



(s) 

depend on the vector and scalar particle number densities Ui and n] , re- 
spectively, that are defined as 



Hi 



dn 

dfli 



is) 



rrii - Si 



fi 



(2^)' - (m, - ^,)2 

with the Fermi-Dirac (Bose-Einstein) distribution function 

Ei — jli 



(51) 
(52) 



m 



exp 



T 



± 1 



or the Maxwell-Boltzmann distribution function 

Ei — jli 



fi{k) = exp 



T 



(53) 



(54) 



depending on the particle statistics. 

The field equations for the mesons are found from the functional fl with 
the help of the Euler-Lagrange equations. They have the form 



mluj 




i 


-5.n)(l 




mlp 


= TpTlp 


i 


-5^n)(l 




mla 










ml6 


= ^sns 









dABi 

duj 
dABi 
dp 



(55) 

(56) 

(57) 
(58) 
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with additional contributions to the source terms due to the dependence of 
the binding energies Bi on the vector meson fields. Finally, the entropy is 
obtained as 



S - ^E-; (2vr)3 



[f,\nf,±ilT f^)\n {It fi)] 



dABi 
dT 



E 



9i dT 



(59) 



with the usual single-particle contribution, a term due to the temperature de- 
pendence of the cluster binding energy shifts ABi and a contribution caused 
by internal excitations of the clusters, i.e. the introduction of the tempera- 
ture dependence of the degeneracy factors, described in Sect. HI In case of 
the Maxwell-Boltzmann statistics for a particle z, the integrand contains only 
the term /j In /j. 

3.2. Scheme of the fugacity expansion 
Introducing the modified fugacity 



Zi 



exp 



mi + Si - Vi 



T 



exp 



S^-V^ 

T 



(60) 



we write the distribution function as 

f^ 



Zi ^ exp 



T 



± 1 



1 ± Zisxp {-y) 



(61) 



with the kinetic energy ei{k) = ^/k'^ + {nii — SiY — {mi — Si) > 0. For <^ 1 
the distribution function can be expanded in a power series in Zi exp (-ej/T). 
The appearing momentum space integrals can be evaluated explicitly |37[ 
with the integral representation of the modified Bessel functions Kjj{x) |38|. 
For the contributions fl38l) of the individual particles to the grand canonical 
potential we find 



A? 



gi f mi - Si 



nrii 



X 



oo 

E 

n=0 



(Tl)^ 



{n + 1)5/2 



3/2 



(n+1) 



(62) 



rrii 



S, 



T 



exp 



(n + 1) 



S^-V^ 

T 
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with the non-relativistic fugacity Zi and the functions 



k,,(x) 



'2x 



TX 



ex:p{x)Ky{x) . 



(63) 



Similarly, the vector density is obtained as 

C \ 3/2 

rrii - Si^ 



9i_ 
A? 



X 



n=0 



rrii 

(n + 1)3/2' 



(64) 



(n + l) 



rrii — Sj 



T 



exp 



(n + 1 



' T 



n+1 



(s) 

For the scalar density n] , k2 has to be replaced with ki. In case of Maxwell- 
Boltzmann statistics, only the n = term remains. Without interaction, i.e. 
Si = Vi = 0, the results for a relativistic Fermi or Bose gas are recovered. 
For T — )• 0, i.e. x = {n + l)(mj — Si)/T — oo, one finds from the asymptotic 
expansion 

ku{x) — 
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1 + 



(^-l)(/x-9) ^ (;,-l)(^-9)(^-25) ^ 



5X 



2!(8x)^ 



3!(8x)3 



(65) 



with /i = 4i/^ the common expressions 



ni 



m+l 



« m=0 



and 



rii 



9i_ 

A3 ^ 

« m=0 



(j + 1)^/^ 



m+l 



(m + 1)3/2' 



(66) 



(67) 



for non-relativistic particles. In this limit, n[^^ = rii, because the scalar den- 
sity differs from the baryon density only in the relativistic description. Since 
the scalar and vector self-energies are themselves functions of the densities, 
a second expansion of the series (IS^ and flMl) is required. 

3.3. Fugacity expansion of the grand canonical potential up to second order 

In order to compare the true series expansion of the grand canonical 
potential Q of the gRMF model with the form f|T8|) in the VEoS approach at 
low densities, we expand f lS7|) up to second order in the fugacities of neutrons 



20 



and protons. In the following we will only consider neutrons, protons and 
deuterons in the density expansion. Nuclei with mass numbers A> 3 do not 
contribute in the second order of expansion since we only consider nucleons 
as basic constituents. The explicit contribution of the deuteron ground state 
appears in the two-nucleon correlation term with 



For sufficiently low nucleon densities, the self-energies of the nucleons are 
approximately linear in the nucleon densities. Hence, the contribution of the 
individual nucleon can be approximated as 




(68) 




k'^ (f ) 3 T 



T 



Zi (69) 



Applying the recursion relation 




(70) 



the expression (l69l) for the nucleons reduces to 



VT 



( 




) 



(71) 



with the abbreviations 




(72) 



(73) 



For the deuteron contribution we have 




(74) 
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without self-energy terms which contribute only in higher order of the fu- 
gacities. Contrary the nucleon self-energies contribute already in the lowest 
order of the fugacities 



Vn 

Sri 



CaiVn + Vp) + CsiVn - Vp), 
CaiVn + Vp) - CsiVn - Vp) 



with coefficients 



Cry 



r^(o) 



(75) 
(76) 
(77) 
(78) 

(79) 



for m = co,p,a,6 that depend on the density dependent meson couplings 
Tm{Q) taken in the limit ^ — )■ 0. The rearrangement term in the vector self- 
energy does not contribute at this level because it is at least quadratic in the 
densities. With the help of the field equations fl55|) -f l58|) . the mesonic contri- 
butions in f l37|l can be expressed as quadratic forms of the nucleon densities. 
Finally, after performing all necessary expansions up to second order in the 
nucleon fugacities, we obtain the resulting grand canonical potential 



fi(T, V,/i„,/ip) 



-TV— 
A3 

-TV— 
A3 



1 A- Z"?""" 

25/2^^2 [^j. 



(80) 



2 A6 



T 



2 ^ 



vg^ 

2 A^ 

p 



frrip 



{C. + Cs) kl 



V T 



(a + Cs) kl 



mp 
T 



+V 



9n9p 

A^A3 



(C. - Cp) k2 



-iC.-Cs)k^{^) k, 



rrir. 



rrij 
~T 



with contributions from free nucleons and their correlation due to statistics, 
the deuteron bound state and the interaction. 
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3.4- Comparison of fugacity expansions 

A comparison of Eq. (ITS]) with the corresponding expansion 
to extract the virial coefficients 6„, bp, b 
The ffist-order coefficients receive a relativistic correction with 

gnk2 



permits 



„„, bpp and bnp of the gRMF model 



bn 
b„ 



T 

rup 

T 



(81) 
(82) 



They depend on the temperature now. In the non-relativistic hmit T /mi — )■ 
the correction will vanish due to the asymptotic expansion of the function 
cf. Eq. ([nS]). Below we will show that these corrections are important 
and should be taken into account for an improved description of low-density 
nuclear matter. Similarly, there is a relativistic modification to the statistical 
corrections in and bpp, i.e. Eqs. ( !T9|) and (120|) . respectively, and to the 
deuteron contribution to (122]) . From the comparison of ( I18p and ( l8Up . three 
independent relations for the channels nn, pp and np remain 



T 



T 



b + 

bpp ~\~ 

2A6 . 

p 



9n 



:ko 2 



T 



25/2 

(C. + Cp) 

25/2 ^ V T 



-kn 



^2 



•nir. 



T 



-^k2 

pp 



2m^ 
2mr 



T 



Einn) 7-(n?i) 
9 1 I I 

EApp) App) 



(83) 



ic. + Cp) ki [^) - (a + cs) ki {^) 



-T 



-"np 



3/2,3/2 



- —k-), I — I exp I — 



An A 

T 



X3 



V T / 



5. 



T 



T 

(npl) j(npl) 



(85) 



(C. - C,) k2 (^) 



(a - Cs) k, 



T 

ki 



91 



nip 
T 



+ 9l 



(npO) j(npQi) 



nip 
T 



that connect the virial integrals with the strengths Cm of the zero-density 
meson-nucleon couplings rm(0) of the gRMF model. The expected correction 
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factors k2{mi/T) due to the relativistic effects were added to the VEoS part 
(expressed through virial integrals I^-') of relations ( 15^ to 



3. 5. Temperature independent limit of consistency conditions 

The consistency relations ( |83|) to (!85|) were derived from a comparison of 
two fugacity expansions that are valid at low densities and finite tempera- 
tures, i.e. in the virial limit. These conditions arise from the contributions 
in the fugacity expansion that are quadratic in the fugacities. We perform 
an expansion of Eqs. fl83l) to fl85l) in powers of small T/rui and keep the 
lowest-order terms which are independent of T. This is equivalent of taking 
the limit T — )■ in fl83|) to fl85l) . In this limit there are no relativistic correc- 
tions, only s-wave contributions are relevant and the limit fl5^ of the integral 
with the scattering lengths can be used. In the gRMF model the nucleon- 
nucleon interaction is identical in the nn and pp systems. Thus, the first 
two equations, (l83l) and (l84l) cannot be considered independent but should 
be combined. Denoting the scattering lengths in the s-wave channels explic- 



^(nn) (pp) 
i5o ' 



itly with al^/, ai'-, a^^^ 



9n = Qp = 

conditions 



2, ^^Z^ 



and a^Pg^ and considering the degeneracy factors 



Avp) 
9^ So 



9^ So 



1 and g^i^J'^ 



3, the two consistency 



Cuj — Ccr 



TT 



a; 



(nn) 



So 



rur, 



+ 



a I 



(pp) 



So 



nir, 



+ 



(np) n„{np) 



16) 



v 



and 



Cp — Cs 




(nn) 
^^So 
m„ 



+ 



App) 
nir, 



(np) q„(np) 



'Si 



(87) 



for the isovector and isoscalar meson couplings are obtained. The left side 
of Eqs. (186|1 to (!87|) corresponds to the gRMF expansion while the right side 
represents the virial expansion in the effective-range approximation. 

With the coupling constants at zero density Tm{0) and the meson masses 
mm the coupling coefficients Cm can be calculated for different RMF models. 
In Table [1] the results are shown for t he g RMF model with parametrization 
DD2 without 5 meson following Ref. pSy and for a new density dependent 
RMF parametrization DD-ME5 from [39| based on ab-initio calculations in 
nuclear matter including a Lorentz scalar isovector 6 meson. 
In the DD2 parametrization the differences of the coupling coefficients are 
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Table 1: Coupling coefficients ((7^ in the DD2 parametrization of the gRMF model [2^ 
and the DD-ME(5 parametrization of Ref. [soj . 



Model 


DD2 




DD-ME5 


meson m 


Cm [fm^] 


Cm [fm^] 


u 


17.25044 


3 


14.539639 


a 


22.639364 


19.443381 


P 


2.483932 


6.017331 


5 






3.462716 



— Ccr = —5.39 fm^ and Cp — Cs = 2.48 fm"^, respectively. In contrast, 



the calculation with the scattering lengths from Ref. [35[ gives C^ — C, 



— 14.15 fm^ and Cp — Cs = —9.61 fm^. It is obvious that the coupling 
coefficients of the gRMF model do not obey the consistency relations to 
guarantee the agreement with the VEoS at low densities. 

For the recently developed parametrization DD-ME5 [39i], the differences 
of the coupling coefficients turn out to be C^ — C^^ = —4.90 fm^ and Cp — Cs = 
2.55 fm^. These numbers are close to those of the DD2 parametrization. 
There are still differences to the values required by the consistency conditions 
( 186|) and (!87|l and a more drastic modification of the low-density couplings 
seems to be necessary. 

In the DD2 parameter fit of the gRMF approach a particular form of 
the density dependence of the coupling functions Tmio) was assumed. The 
choice was motivated by results from Dirac-Brueckner calculations of nuclear 
matter [iif. However, the values at zero density are found from a not well 
constrained extrapolation to small densities. Only near the saturation density 
of symmetric nuclear matter they are well determined from model fits to the 
properties of atomic nuclei. Thus, in general, one could modify the density 
dependence of Tmio) for all mesons m at low densities in order to satisfy the 
constraints (186!) and ( 187|) . Since Cp — Cs has to be negative, it is clear that 
the standard u, a and p mesons are not sufficient and a 6 meson is needed 
in this case. 

Since Eqs. (ISBj) and (157|) give only two relations for four coupling co- 
efficients Cm, obtained from the consistency with the virial expansion, two 
further constraints are required to fix all couplings unambiguously. Hence, 
in addition to the virial constraints, we tried to fit the vacuum couplings 
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coefficients T^iO) by describing low-energy nucleon-nucleon scattering, i.e. 
experimentally determined NN phase shifts [ssj, simultaneously satisfying 
equations f l86|) and f l87|) . The corresponding Schrodinger equation in coordi- 
nate space was solved using a, u, p and 6 mesons with a one-boson exchange 
potential 40| in order to fit the coupling constants of the nucleon-meson 
interaction and to describe the exp erimentally known phase shifts and effec- 



tive range parameters given in [35[, |4l|. The gradient terms in the potential 



lead to unphysical singularities at r — )■ 0, thus to avoid these singularities a 



certain cutoff was introduced as in j42|. Coupling constants derived in this 



way would serve as an additional constraint at zero densities, improving the 
density dependence of the coupling functions Tmig)- However, even when the 
contribution of the pion was included, that does not appear in the mean-field 
approximation, it was impossible to find a reasonable parametrization. Quite 
drastic modifications of the low-density meson-nucleon couplings and meson 
masses were required. This failure is not very surprising because the number 
of meson fields in the RMF model is much smaller than that in realistic NN 
potentials resulting in a smaller number of degrees of freedom. In addition, 
the meson fields in the RMF model do not represent real physical mesons 
and their masses, that are used as adjustable parameters, e.g. in case of the 
a meson, are not identical to the true meson masses. Hence, we do not follow 
this approach but extend the gRMF model further. 



4. Extension of the gRMF approach 

The failure of the gRMF model to describe the virial limit suggests that 
a further extension of the approach is inevitable. In addition to the bound 
state two-body correlation, the deuteron, we introduce additional degrees of 
freedom in the grand canonical potential that represent the possible two- 
nucleon states, nn, pp and np, in the continuum. In the following, we will 
consider the isospin and 1 channels for the np states separately. Each of 
the two-body channels is represented by a temperature dependent resonance 
energy {Enn, Epp, Enpo, Enpi) and a corresponding effective degeneracy factor 
{9nn, Qppi QnpOi Qnpi)- In general, the resonance energies and degeneracy 
factors do not necessarily have to be the same as the resonance energies of 
Eq. or the standard degeneracy factors in Eqs. (IT^ to Two-body 
correlations in the continuum are considered as two-body clusters with mass 

rhij = NijUin + ZijUip + Eij (88) 
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with a temperature dependent resonance energy Eij. We assume that the 
energy shift of a cluster ABij that appears in the scalar self-energy (1441) 
is identical to the binding energy shift of the deuteron AB^ for all contin- 
uum channels. Of course, this most simple choice can be substituted by 
improved descriptions in the future. The specific functional form of the 
medium-dependent energy shift will affect the form of the transition from 
very low densities, where the VEoS is applicable, to the high-density regime. 
However, for the low-density matching of the VEoS and gRMF model it is 
irrelevant. 

The total rest mass of the cluster in the medium is given by 



Mi, 



m 



1-3 



+ AB, 



(89) 



In case of the virial approach, the appearing masses are given by iriij = 
NijUin + Zijirip + Eij and Mij = rriij, respectively. We include the effect of 
the resonance energies in the relativistic correction through the k2 function 
and the definition of the thermal wave lengths Xij = [2%/ {MijT)^/^ and sim- 
ilarly for Xij. With these modifications we can rederive consistency relations 
between this extended gRMF model and the VEoS. 

The virial integrals fl23|) in the virial coefficients are represented by an 
exponential of the effective resonance energy as defined in (^^. The two 
np channels are separated by requiring for the np isospin 1 channel a similar 
form as for the nn and pp channels. Comparing the fugacity expansion of the 
gRMF model with 2-body correlations with the VEoS at low densities, we 
finally obtain an extension of the relations (183|) to ( |85|) , where an additional 
term with temperature dependent resonance energy Eij{T) appears in the 
right-hand side of the Eqs. (MI]) to ( IHT]) . i.e. 




(90) 
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for the channels ij = nn, pp, npl and 



T 



'npO 




(91) 



for the channel npO. These consistency relations leave sufficient freedom 
to reproduce the exact VEoS at low densities in the gRMF model. There 
are different possibilities to achieve this goal depending on the choice of the 
zero-density couplings and on the choice of the effective resonance energies 
and degeneracy factors that represent the continuum contributions in the 
extended gRMF model. The structure on the left-hand side shares some 
similarities with the generalized Beth-Uhlenbeck approach in Section 12.41 
One part of the effects of the two-body correlations is included in the explicit 
continuum contributions and the remaining part is covered by the mean-field 
terms. There is the freedom to shift the correlation strength between these 
terms rather arbitrarily. In the following, a few choices will be discussed. 

In Section [3751 the case without introducing explicit continuum channels in 
the gRMF model, i.e. shifting the full correlation strength into a redefinition 
of the vacuum couplings, turned out to be impossible. The other extreme 
is to represent the correlations by the introduction of the effective contin- 
uum states described by temperature dependent resonance energies Eij with 
medium corrections ABij as described above. In this case the contribution of 
the virial integral in the left-hand side of the Eqs. ( M?]) and (IHT]) is fully rep- 
resented by the the resonance energies Eij in the first term in the right-hand 
side of these equations. Therefore the mean-field part has to vanish in the 
non-relativistic limit, leading to C^j — = and Cp — Cs = 0. The coupling 
coefficients Cm in relations f l86|) to f l87|) were expressed through the virial in- 
tegral taken in the effective range approximation. In the new relations, given 
above, the virial integral and thus the scattering lengths are implemented 
in the definition of the resonance energies Eij{T). Therefore, neglecting the 

term proportional to g^^^^ exp (—Eij/Tj in the relations f l90|) and fl9T]) . we 
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recover equations (1551) to (153]) . leading to the consistency conditions (15^ to 

dHZD- 

In order to fix the four couplings Cm unambiguously, we again tried to 
describe low-energy nucleon-nucleon scattering, by solving the Schrodinger 
equation with one-boson exchange potentials as described in Section 13.5! 
simultaneously with the two conditions on the couplings C^^ — = 0, 
Cp — Cs = 0. We observed that the required coupling constants at zero 
density are drastically different in comparison with the previous values of 
the parametrization DD2 or DD-ME(5 leading to problems such as a rather 
unphysical behavior of the EoS at low densities by similar reasons as in the 
first attempt at the end of Section 13.51 Consequently, this approach was dis- 
carded. It was not possible to describe both the NN scattering data and the 
VEoS simultaneously in this second choice for representing the continuum 
correlations. 

Therefore the most straightforward way to fulfill the consistency relations 
f l90|) and fl9T]) is to leave the density dependence of meson-nucleon couplings 
untouched as given by the gRMF model parametrization and to turn down 
the idea of an additional fit to NN scattering data, thus concentrating only 
on the reproduction of the VEoS. The effects of the two-body correlations are 
incorporated in the effective resonance energies that are chosen as identical 
to those of the VEoS, i.e. Eij = Eij. Then, the consistency relations ( 190|) and 
(PTj) are seen as the defining equations for the effective degeneracy factors 
g^j^^ that now depend on temperature. In the non-relativistic limit, these 
equations simplify as 
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(92) 
(93) 
(94) 
(95) 



where 



Cl — C^ — Ccr + Cp — Cs , 

Co = C^ — Co- ~ S {Cp — Cs 



(96) 
(97) 
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Figure 3: Temperature dependence of the effective degeneracy factors .gi°f •* and g)^^f^ with, 
and without relativistic corrections. 



,(off) 



are the relevant couplings in the two isospin channels. If Ci = Co = as was 
discussed before as a possible condition, we recover the standard degeneracy 
factors Qij in the s-wave. The degeneracy factors on the right-hand side are 
gnn = Qpp = gnpi = 1 and gnpo = -3, respectively. The negative sign for (jnpo 
is related to the positive s-wave scattering length in this channel. Note that 
in the NSE description, cf. Subsection l2.3[ the temperature dependence of the 
degeneracy factors follows from the inclusion of excited states of the nuclei. 
In our approach, it represents the contribution of the continuum states and 
thus is very similar. 

In Fig.[3]the temperature dependence of the effective degeneracy factors is 
depicted for the scattering correlations in the nn and npO channels. Because 
the curves for the pp and npl channels are practically identical to the nn 
scattering case, they are not shown. The fully relativistic results as they 
follow from Eqs. ( 190|) and (!9T]) are compared with their non-relativistic limits 
( 192|) and ( 195|) . respectively. We see that the relativistic corrections caused 
by the ki and /c2 functions are minor at low temperatures but grow with 
the increase of temperature as expected, being larger for the case of nn 
scattering. The large relativistic correction in the effective degeneracy factors 
is caused by an accidental change of sign in the factor containing the coupling 
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coefficients in Eq. (M?]) . 

In the zero-temperature limit, the effective degeneracy factors are given 

by 

9fm = 9, - , (98) 

for the channels ij = nn,pp, npl and npO, respectively, when the approxima- 
tion (l32l) with the scattering length is used. Thus, there is still a difference 
as compared to gij = 1 and c/npo = —3, respectively, when the correction due 
to the meson couplings is not taken into account. This is a clear indication 
that the mean field effectively takes over a part of the correlations. 

5. Neutron matter 

In nuclear matter with arbitrary neutron to proton ratio, the two-body 
correlations in all nucleon-nucleon channels contribute to the thermodynamic 
properties and, in particular, the deuteron bound state dominates at low tem- 
peratures. The case is different for pure neutron matter where only the nn 
channel without a bound state is relevant. Thus, neutron matter is ideally 
suited to demonstrate the effects of the continuum correlations in the ex- 
tension of the gRMF model. Before the general case of finite temperatures 
is discussed, we first consider in brief the case T = in the limit of small 
densities. 

5.1. Zero temperature limit 

At zero temperature the Fermi momentum kp^ defined through the neu- 
tron density n„ = /c|^^/(37r^) sets the scale for all results. Instead of an 
expansion in powers of the neutron fugacity at finite temperature, a series 
expansion in powers of kp^ is the relevant method to study the low-density 



behavior of the EoS. In Ref. [43| it was shown that the energy per neutron 



E/N (without rest mass) at very small kp^ can be expanded in a power series 



f = Ei,U (100) 



with the power series 



10 4 

e = l + ^C + ^(ll-21n2)C^ + ... (101) 
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in the dimensionless parameter ( = ai^^Jkp^. Here i^free = 3fc|.^/(10m„) is the 
energy per neutron of a non-interacting Fermi gas that defines the relevant 
energy scale. Because of the unnaturally large negative nn scattering length 
a^^^J !^ —18.8 fm js^, the radius of convergence is unfortunately very small, 

i.e. rin < 2437r/|10a5"^|3 ^ 1.1 ■ 10"^ fm'^. 

In the DD-RMF approach without effective contributions of two-body 
correlations, an expansion of the energy per neutron in powers of the Fermi 



momentum (see [Appendix B[) leads to 



E 

= -f'free^RMF (102) 

with 

^RMF = 1 + TT^iC^-C^ + Cp- Cs) m^kF„ + . . . . (103) 
A comparison with Eq. f llOip gives the relation 

C^-C^ + Cp-Cs = ^al™) . (104) 

This condition coincides with condition ( |83|) for neutron matter, although it is 
derived in the low kp^ / nin limit for T = whereas condition (I83p in the T — >■ 
limit of the finite-temperature virial expansion with vanishing convergence 
radius in density. From a physical point of view, this coincidence is not 
surprising since both approaches only use two-body scattering information. 
However, it is gratifying to find the coincidence from the two very different 
approaches. We find the values of —2.91 fm^ and —2.35 fm^ for — + 
Cp — Cs in the DD2 and DD-ME5 parametrizations, respectively. However, 
these are much smaller in modulus as compared to the required value of 
27rai""V^n = -24.83 fm^. 



A particular situation arises in the unitary limit 44 1, i.e. if a^^^J ap- 



proaches — oo. Then the series expansion f llOll) can no longer be applied 
since the radius of convergence shrinks to zero. In fact, the energy per nu- 
cleon should scale as the energy per neutron of a non-interacting Fermi gas 
i?free with a Universal constant ^ independent of kp^. The parameter ^rmf is 
independent of kp^ in first-order only if the combination C^^ — C^^ + Cp — Cs 
diverges as kp^. Hence a particular density dependence of at least one meson- 
nucleon coupling is required in this case. 
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In the nonlinear RMF model of Ref. [12j with nonlinear self-interactions 
of the scalar meson, a density dependence of the o meson coupling was intro- 
duced for densities lower than a particularly chosen transition density with 
the aim to reproduce the energy per neutron for unitary neutron matter as- 
suming ^ = 0.44 45|. This condition required a divergence of V f,{n) oc n~^l'° 



that is consistent with the expectation from Eq. fll03p . However, mixing the 
nonlinear RMF approach and the RMF model with density dependent and 
divergent couplings does not seem to be very natural. In our approach, we 
do not aim to describe the neutron matter EoS at zero temperature and low 
densities as a unitary Fermi gas (FG) but require consistency with the VEoS 
for T > 0. 

5.2. Finite temperatures 

We choose T = 4 MeV and T = 10 MeV as representative examples which 
illustrate the main effects due to interactions and correlations. Different 
approximations in the description will be compared, ranging from the ideal 
Fermi gas via the virial approach to the gRMF model. Nucleons are treated as 
Fermions and clusters as Maxwell-Boltzmann particles because they appear 
at high temperatures with low abundancies when deviations from the correct 
Bose-Einstein statistics are negligible. 

There are two quantities that exhibit finite low-density limits and hence 
are very advantageous for a comparison of the models: the ratio of the pres- 
sure over total particle number density p/n and of the internal energy per 
baryon E/N = e/n — rUn (without the rest mass of the neutron) with the 
energy density 

e = ^{Q + TS)+ Mi (105) 

i=n,nn 

that can be calculated from the grand canonical potential Q, the entropy 
S and the relativistic chemical potentials jli. The sum contains the contri- 
bution of the neutrons and that of the correlated two-neutron continuum. 
Correspondingly, the total particle number density n is the sum 

n = n„ + 2nnn (106) 

and the mass fraction of correlated two-neutron states is defined by 
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Figure 4: Ratio of pressure over total particle number density, p/n, of neutron matter as 
a function of the total density n for temperatures of T = 4 MeV (left) and T = 10 MeV 
(right). Vertical dotted lines indicate the density where nA^ = 1/10. 



For an ideal gas with Maxwell-Boltzmann statistics and non-relativistic kine- 
matics, we have the simple result that p/n = T and E/N = 3T/2 are inde- 
pendent of the density of the system and trivially X„„ = 0. Deviations from 
these values indicate the effects of correlations and interactions. 

5.2.1. Low densities 

In Figures m and [5] the two quantities p/n and E/N, respectively, are 
depicted for the two selected temperatures as a function of the total particle 
number density n in different theoretical approaches. In case of the rela- 
tivistic Fermi gas (see FG-curve), effects of the Pauli principle are taken into 
account leading to an increase of the pressure and of the energy per neutron 
as compared to the ideal Boltzmann gas. The limit lim„_j,o(p/^) = T is not 
affected by statistical corrections or relativistic kinematics since the k2 fac- 
tors cancel in the lowest order of the fugacity expansion. This is easily seen 
considering the ratio p = —Vt/iVn) using Eqs. (1621) and (1641) . In contrast, the 
relativistic correction factor k2 in Eq. (IMl) modifies the relation between the 
neutron density and neutron chemical potential appearing in Eq. (11051) and 
thus lim„^o(-£'/^) > 3T/2 in figure [51 The relativistic corrections become 
larger with increasing temperature. The shift of E/N ai zero density can be 
estimated as 2T^/m„. 
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Figure 5: Internal energy per baryon (without contribution of the neutron rest mass) in 
neutron matter as a function of the total particle number density n for temperatures of 
T = 4 MeV (left) and T = 10 MeV (right). Vertical dotted lines indicate the density 
where rtAfj = 1/10. 



The VEoS predicts a dependence oi p/n and E/N on the density with 
a negative slope. This is the effect of the correlations induced by the nn 
interaction. The vertical lines in Figures H] and |5] denote the density n where 
n\\ = 1/10. At densities above this value, higher-order contributions to 
the VEoS, that are not considered in the fugacity expansion up to second 
order, can be expected to contribute significantly. Figure [5] also demonstrates 
that the virial corrections at low densities are smaller than the relativistic 
correction that leads to an overall shift. 

The curves of the DD-RMF model with DD2 parametrization without 
correlations lie between the VEoS and the FG results. They do not show 
the correct dependence given by the VEoS at low densities. When the nn 
correlations are taken into account in the gRMF model with the quadratic 
form of the energy shift (see next subsection) the low-density behavior of the 
VEoS is nicely reproduced. Only at higher densities, medium effects, that 
are not incorporated in the VEoS, lead to a deviation. The precise density 
dependence of the deviation will depend on the choice of the functional form 
of the energy shift, but the agreement in the low-density limit is not affected. 
Obviously, deviations of the gRMF EoS from the VEoS start to become more 
important with increasing density at lower temperatures. 
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Figure 6: Ratio of pressure over total particle number density, p/n, of neutron matter as 
a function of the density n for temperatures of T = 4 MeV (left) and T = 10 MeV (right) 
in different models. See text for details. 



It is also worthwhile to compare the predictions of the VEOS, the original 
DD-RMF and gRMF models in the DD2 parametrization with the results 
of other approaches used in astrophysical applications. We examine two 
other RMF models that employ nonlinear self-interactions of the mesons: 
the model of Shen, Horowitz and O'Connor (SHO) [l5[ with the FSUG- 
old parametrization and the model of Shen, Toki, Oyamatsu and Sumiyoshi 
(STOS) del with the parameter set TMl. In addition, the non-relativistic 
Lattimer-Swesty EoS with incompressibility K = 220 MeV (LS 220) is con- 
sidered. 

In Figure [6] the density dependence of the quantity p/n at T = 4 MeV 
and T = 10 MeV is depicted for all models below 0.001 fm~'^. The curves of 
the STOS model, that does not include two-body nn correlations, are sur- 
prisingly close to the exact VEoS line for T = 10 MeV, but do not reproduce 
it exactly. The deviations are larger at T = 4 MeV and the close agreement 
seems to be accidental for T = 10 MeV. The SHO approach with the FSUG- 
old parametrization claims to be constructed such that the case of unitary 
neutron matter is reproduced at low densities by introducing a particular 
density dependent coupling of the cr meson, see the discussion at the end of 
Subsection 15.11 A large deviation from the VEoS result can be seen in Fig. 
Eland the correct low-density limit for p/n is not reproduced. Furthermore, 
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Figure 7: Ratio of pressure over total particle number density, p/n, of symmetric nuclear 
matter as a function of the total density n for temperatures of T = 4 MeV (left) and 
T = 10 MeV (right). 



the tabulated data exhibit some oscillations that are probably related to the 



choice of the interpolation procedure j46|. The LS EoS shows a much larger 



negative slope of p/n as a function of n as compared to the other models but 
reproduces the correct ideal gas limit. 

6. Symmetric nuclear matter 

In contrast to neutron matter, where only neutron scattering correlations 
contribute to the thermodjTiamical quantities, all correlations of neutrons 
and protons in scattering and bound states should be taken into account 
in symmetric matter. At low densities, two-body correlations will be most 
important, but with increasing density also many-body correlations, in par- 
ticular the appearance of clusters, i.e. many-body bound states, are rele- 
vant. Presently we include contributions from nni^So), np{^So), pp{^So) and 
np{^Si) scattering channels and clusters with A < 4 in our calculations. In 
the gRMF model light clusters (deuteron, triton, helion and a-particle) are 
introduced as additional degrees of freedom with temperature and density 



dependent binding energy shifts given in Ref. j25|. Because we do not con 



sider the formation of nuclei with mass numbers A > 4 the present model can 
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Figure 8: Internal energy per baryon (without contribution of the neutron rest mass) 
in symmetric nuclear matter as a function of the total particle number density n for 
temperatures of T = 4 MeV (left) and T = 10 MeV (right). 



only be applied to rather low densities where the fraction of heavier clusters 
can be neglected. 

6.1. Low densities 

Similar as in the case of neutron matter, Figures [7] and [8] depict the quan- 
tities p/n and E/N, calculated in different approaches, for T = 4 MeV and 
T = 10 MeV as a function of the total particle number density n. The 
extended gRMF model very well reproduces the VEoS at low densities, devi- 
ating from it with increasing density and lower temperature due to medium 
effects. Note that in the depicted results of the VEoS calculation the same 
nuclei are considered as in the gRMF model. It is also worth noticing that the 
standard RMF calculation without clusters and the extended gRMF model 
differ substantially. Again, we see the effect of the relativistic corrections on 
the internal energy per baryon E/N in comparison with the ideal gas limit. 
In contrast to the neutron matter case, where scattering correlations are es- 
sential for reproducing the VEoS, in symmetric matter the main contribution 
is caused by the appearance of bound states with positive binding energy. 
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Figure 9: Particle fractions Xi of nucleons and light clusters in symmetric nuclear matter 
as a function of the density n for temperatures of T = 4 MeV (left) and T = 10 MeV 
(right). Solid lines correspond to the VEoS and dotted lines to the gRMF model. 



6.2. Composition 

It is worthwhile to study the detailed composition of symmetric nuclear 
matter, comparing the fractions of various particles. In Figure [9] particle 
fractions Xi = AiUi/n of nucleons and light clusters are shown for T = 4 
MeV and T = 10 MeV for gRMF (dotted lines) and VEoS (solid lines) 
models. The deuteron fraction contains the contribution of the bound 
state and the isospin singlet np scattering channel. As is seen from the 
figures, in symmetric matter higher mass clusters become more important 
with increasing densities. One notes that the particle fraction for deuterons 
in the VEoS and gRMF models are substantially different, especially in case 
for T = 10 MeV. This behavior is caused by the fact that in the gRMF model 
part of the two-body correlation strength is shifted to the mean field, thus 
reducing the contribution of the original two-body correlations. The same 
effect will be observed in case of neutron matter for the particular case of nn 
correlations, see Section [71 We also notice different slopes of the curves for 
light clusters, with the a particle having the steepest inclination. The reason 
is that at low total densities cluster densities are proportional to fugacities 
in a certain power, e.g. ~ z^, ~ and ~ with z = Zn ^ Zp in 
symmetric nuclear matter. The overall scaling of the cluster fractions at low 
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densities is determined by the cluster binding energies. At densities n above 
10"^ fm-3 and lO'^ fni'^ for T = 4 MeV and T = 10 MeV, respectively, 
heavier nuclei will contribute significantly to the composition in symmetric 
nuclear matter and they have to be incorporated in the model calculations. 
Hence, it is not reasonable to discuss the transition to higher densities in the 
present description of symmetric nuclear matter. The case of pure neutron 
matter is presented in the next section. 



7. Higher densities 

In subsection 15.2.11 it was demonstrated that the gRMF approach with 
effective nn scattering correlations reproduces the low-density limit of the 
VEoS for thermodynamical quantities such asp/n and E/N. With increasing 
density, the VEoS approach is no longer applicable and a smooth transition 
of the gRMF predictions to the DD-RMF results with neutrons as quasi- 
particles is expected. The details of this transition are affected by several 
ingredients of the gRMF model that are not constrained by the low-density 
expansion: the strength of the cluster- meson couplings (H5|) . the shift ABnn 
of the effective resonance energy E^n and possible contributions from three, 
four, . . . and many-neutron correlations. In the following, only the varia- 
tion of the transition with different choices of ABnn will be discussed. The 
strength of the cluster-meson couplings is kept as given in the original gRMF 
model and described in Section |3l 

When the total density of neutron matter increases, the contribution 
of the nn cluster in the gRMF model will be affected by the applied shift 
to the resonance energy that represents the continuum correlations. The 
functional dependence of the energy shift ABi of a cluster i in (jH]) with 
temperature T and the meson fields u and a is not determined by the low- 
density considerations. Here, we explore three different choices. Similar as 



in Ref. 25 



we write 

AB, = f[nf^]SBi{T) (108) 
with a function / that depends on the effective density 

(eff) ml Nj-Z, ml 

and a temperature dependent factor 5B.i{T). At low effective densities, ABi 
should be linear in nf^^ and an obvious choice is / = nf^\ In Ref. 25| the 
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Figure 10: Mass fraction X„„ of the two-neutron correlation state in neutron matter as a 
function of the density n for temperatures of T = 4 MeV (left) and T = 10 MeV (right). 



quadratic form 



/ 



n, 



(cff) 



1 n 
1 + - 



(eff) 



9 (0) 



(110) 



was used for the light clusters in order to obtain a stronger suppression of 
the cluster abundancies with increasing density. The density scale nf' = 



b'i"'^^/5B,{T) is set by the vacuum binding energy Bl . Another possible 
choice is the pole form 



(vac) 



/ 



(eflf) 

^sat 



^sat 



n. 



(eflf) 



;iii) 



for nf'^'^ smaller than the saturation density rigat of the gRMF model resulting 
in a complete dissolution of the cluster when nsat is approached from below 
because lim„_j.„^^^ f = oo. 

The evolution of the two-neutron mass fraction 



Xr, 



2n„„/ (n„ + 2n 

r 



(112) 



with increasing total density is depicted in Figure [TO] for the VEoS and the 
gRMF model with the linear, quadratic and pole form of the energy shift 
ABnn, respectively. We assume 6Bnn(T) = 6Ba(T) and set the density scale 
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Figure 11: Ratio of pressure over total particle number density, p/n, of neutron matter as 
a function of the density n for temperatures of T = 4 MeV (left) and T = 10 MeV (right). 
Vertical dotted lines indicate the density where nA^ — 1. 



Unl = fid = B^J^^' /6Bd(T) with the deuteron values for the quadratic de- 
pendence. The nn fraction in the VEoS model rises monotonously with the 
total density reaching unrealistically high values much beyond the range of 
applicability of the approach. At low densities, the nn fractions for the dif- 
ferent choices of the energy shift in the gRMF model agree perfectly with 
each other. They exhibit the same slope as the VEoS result. In general, the 
nn mass fraction at low densities is larger for the lower temperature but the 
two-neutron cluster dissolves earlier with increasing total density. The max- 
imum mass fraction and range of cluster dissolution depends sensitively on 
the form of the energy shift AB„„. There are substantial variations that need 
to be constrained in future investigations. On an absolute scale, the gRMF 
predictions for X„„ are substantially smaller than those of the VEoS at low 
densities. This difference is caused by the fact that in the gRMF approach 
the correlations of quasiparticles are considered and part of the correlation 
strength is contained in the self-energies, cf. the generalized Beth-Uhlenbeck 
approach in Subsection 12. 4[ The distribution of correlations between the ex- 
plicit contribution from the cluster state and the implicit contribution via 
the self-energies depends on the nucleon-meson couplings at zero density of 
the particular gRMF parametrization. 
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Figure 12: Internal energy per baryon (without contribution of the neutron rest mass) 
in neutron matter as a function of the total particle number density n for temperatures 
T = 4 MeV (left) and T = 10 MeV (right). Vertical dotted lines indicate the density 
where nX't — 1. 



The dependence of the quantities p/n and E/N m neutron matter on the 
total density is shown in Figures [H] and [I2] for temperatures T = 4 MeV 
and T = 10 MeV and a wider range of densities compared to that shown in 
Figures S] and [51 The vertical lines in Figures [11] and [12] denote the density 
n where n\\ = 1. In the low-density limit, all gRMF calculations reproduce 
the VEoS predictions by construction but deviate from the DD-RMF model 
that does not take cluster formation into account. At higher densities, the 
VEoS fails to predict the strong increase of the pressure and energy per neu- 
tron caused by the short-range repulsive nn interaction. The transition of 
the gRMF results at low densities to the DD-RMF curve at higher densities 
substantially depends on the choice of the energy shift Ai?„„ for the effective 
resonance energy Enn- A distinctive bump in p/n and E/N appears that is 
correlated with the sudden dissolution of the two-body clusters as depicted 
in Figure [10] This feature was already observed in Ref. |25| for symmetric 
nuclear matter including only the bound states of light clusters. The origin 
of this structure is related to the contribution with the derivatives of the 
energy shifts in the u and p-meson field equations (!55]l and (156]) . The func- 
tional form of the density dependence and absolute scale of the meson-cluster 
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coupling strengths from fl45l) and fl46l) will also have an impact on the 
detailed form of the transition from the low-density limit to higher densities. 
In this work, only the most simple choice of the factor Qim, proportional to 
the number of nucleons in the cluster, was examined. Further investigations 
are needed to fix the cluster-meson couplings and the energy shifts less am- 
biguously. Correlations beyond two-neutron states, that are not considered 
for neutron matter in the present approach, could also modify the features 
in the quantities p/n and E/N. 

8. Conclusions 

In this paper, an extension of the generalized relativistic mean-field model 
with density dependent couplings was developed by requiring the consistency 
of the finite-temperature equation of state at low densities with the virial 
equation of state that is a model independent benchmark in this limit. For 
this purpose, new degrees of freedom were introduced in the generalized RMF 
approach that represent two-nucleon correlations in the continuum. These 
clusters can be considered as quasiparticles like bound nuclei. They are 
characterized by medium-dependent effective resonance energies with tem- 
perature dependent effective degeneracy factors. 

From the comparison of the fugacity expansions of both the relativistic 
mean-field and virial equation of state models, consistency relations were de- 
rived that contain the nucleon-nucleon scattering phase shifts, the relativistic 
mean-field nucleon-meson coupling constants, the resonance energies and ef- 
fective degeneracy factors of the clusters. Various limits were investigated for 
the choice of the relevant parameter functions and their dependence on the 
thermodynamical quantities. For a successful application of the approach, 
the density dependence of the original relativistic mean field model was kept 
and not modified at low densities. The effective resonance energies were taken 
as given by the calculation with the scattering phase shifts and the effective 
degeneracy factors were assumed to be temperature dependent similar as in 
the case of the nuclear statistical equilibrium model with thermally excited 
nuclei. 

The example of neutron matter was studied in detail for different temper- 
atures. Relativistic effects were found to become important with increasing 
temperature even at very low densities. The extended generalized relativistic 
mean-field model smoothly interpolates between the correct low- and high- 
density limits describing the dissolution of the clusters. However, the precise 
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form of the transition depends on the couphng strength of the clusters to the 
meson fields and the energy shift of the resonance energies. These quantities 
are not fixed by the low-density constraints and the consequences of different 
choices were investigated. This point deserves more studies in the future. 
In the case of symmetric nuclear matter the formation of many-body bound 
states is crucial for a realistic description of correlations. 

The proposed extension of the generalized relativistic mean-field model is 
rather general and can also be applied to other mean-field approaches that 
aim at describing the equation of state of nuclear matter. In the present 
work, only two-particle correlations were considered. Bound states of light 



clusters with mass numbers A < 4 can be readily included as in Ref. [25|. A 



further extension to heavier nuclei with medium- dependent binding energies 
and finite-temperature excitations is also possible and corresponding work 
is in progress. For astrophysical applications, further contributions to the 
EoS than those considered in this work have to be included in the model. 
E.g. electrons, muons, photons, mesons etc. need to be considered which also 
change the behavior of various thermodynamical quantities in the low-density 
limit. Finally, equation of state tables with the thermodynamical properties 
and the composition of nuclear matter for a broad range in temperature, 
density and proton-neutron asymmetry can be generated for astrophysical 
simulations of, e.g., core-collapse supernovae. 
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Appendix A. Virial equation in our model and in Ref. 
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The authors of Ref. 18|] use a slightly different definition of some quan- 



tities in their approach to the non-relativistic VEoS compared to the ones 
we use in Section |2J In order to facilitate the comparison, we indicate the 
correspondence of the two formulations. Quantities of Ref. 18(| are indicated 
by a " in the following. 

The single-particle partitions functions are defined in the same way, i.e. 
Qi = Qi for nucleons and a-particles, however, there is a small difference in 
the thermal wavelengths because the neutron and proton mass are assumed 
to be equal 



m 



nir. 



rhp in Ref. 



18 



and the a-particle mass is set to 
rha = 4:m without considering the binding energy as in Eq. f HU]) . Similarly, 
for the non-relativistic chemical potentials the relation jla = '^f^n + '^'l^p is used 
instead of fia = 2/in + 2/ip — Ba in our case. Nevertheless, the fugacities are 
identical. The main differences stem from the definition of the many-body 
partition functions. In Eq. ([T]) we place factors l/n\ in front of the n-body 



terms. Consequently Qij = Qij/2. Comparing Eq. with Eq. (11) in [18 
we identify 



bpn/'^ 



and 



-'ap 



(A.l) 



(A.2) 



In Eq. ffT6l) we use the cm. energy E as the integration variable. In contrast, 
in Ref. [l8| the laboratory energies E = 2E are used in Eqs. (19), (22) and 
(24) and the integrals are transformed with the help of a partial integration 
with respect to E. Noting that 6p„ = 6nuc — in Ref. 31, the formulas given 



there for bpn and ba are consistent with the relations (lA.ip and (lA.2p . For 



the virial coefficient ban, the authors of Ref. [18j use the nucleon laboratory 



18 



energy E = 5E/A as integration variable. The expression (26) in Ref. 
however, is a factor two too large to be consistent with the relation given in 
flA.2|) . This discrepancy was already noted in Ref. 28 . 



Appendix B. Zero temperature low-density limit in the gRMF 
model 

In the case of pure neutron matter at zero temperature, all relevant ther- 
modynamical quantities can be represented analytically as a function of the 
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Fermi momentum kp^. The energy density e of pure neutron matter without 
contributions of the rest mass reads, cf., e.g., Ref. |26| . 



3 / 1 

4 V + ("^" ~ 'S'„)2 n„ + -(m„ - Sn) n^f) 



+ - 



mi 



2 1 



mi 



where the scalar neutron density is given by 

3 



n. 



with the function 



fix) = xVl + — In + vT+~x^ 



(B.l) 



\n 



{s)l2 



(B.2) 



(B.3) 



that depends on the dimensionless parameter x = kp^/{mn — Sn)- We define 
the derivative of the function ( IB. 31) . /' = 2x^/ a/I + x^, and use the expansion 



1 + 2) 



k=0 



r(a + i) 



A;!r(a + 1 - k) 



to rewrite equation f lB.2|) after integration as 

2 



(^) 



1-1 

10 



9 

^ 56 



Zp 



+ 



(B.4) 



(B.5) 



Substituting nii^ in ( IB. 11) by ( IB.SP and doing subsequent expansions in powers 
of kp^ we arrive at expression fllOOp for the energy per neutron. 
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